NASA/TM— 2002-211696 


AIAA-2002-2419 



A Parametric Study of Fine-Scale 
Turbulence Mixing Noise 

Abbas Khavaran 

QSS Group, Inc., Cleveland, Ohio 

James Bridges 

Glenn Research Center, Cleveland, Ohio 
Jonathan B. Freund 

University of Illinois at Urbana-Champaign, Champaign, Illinois 


July 2002 



The NASA STI Program Office ... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key part 
in helping NASA maintain this important role. 

The NASA STI Program Office is operated by 
Langley Research Center, the Lead Center for 
NASA's scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. 
The Program Office is also NASA's institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant 
phase of research that present the results of 
NASA programs and include extensive data 
or theoretical analysis. Includes compilations 
of significant scientific and technical data and 
information deemed to be of continuing 
reference value. NASA's counterpart of peer- 
reviewed formal professional papers but 
has less stringent limitations on manuscript 
length and extent of graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific 
and technical findings that are preliminary or 
of specialized interest, e.g., quick release 
reports, working papers, and bibliographies 
that contain minimal annotation. Does not 
contain extensive analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical 
conferences, symposia, seminars, or other 
meetings sponsored or cosponsored by 
NASA. 

• SPECIAL PUBLICATION. Scientific, 
technical, or historical information from 
NASA programs, projects, and missions, 
often concerned with subjects having 
substantial public interest. 

• TECHNICAL TRANSLATION. English- 
language translations of foreign scientific 
and technical material pertinent to NASA's 
mission. 

Specialized services that complement the STI 
Program Office's diverse offerings include 
creating custom thesauri, building customized 
data bases, organizing and publishing research 
results . . . even providing videos. 

For more information about the NASA STI 
Program Office, see the following: 

• Access the NASA STI Program Home Page 
at h ttp:lhviviu.s ti. nns a. gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access 
Help Desk at 301-621-0134 

• Telephone the NASA Access Help Desk at 
301-621-0390 

• Write to: 

NASA Access Help Desk 

NASA Center for AeroSpace Information 

7121 Standard Drive 

Hanover, MD 21076 



NASA/TM— 2002-211696 


AIAA-2002-2419 



A Parametric Study of Fine-Scale 
Turbulence Mixing Noise 

Abbas Khavaran 

QSS Group, Inc., Cleveland, Ohio 

James Bridges 

Glenn Research Center, Cleveland, Ohio 
Jonathan B. Freund 

University of Illinois at Urbana-Champaign, Champaign, Illinois 


Prepared for the 

Eighth Aeroacoustics Conference 

cosponsored by the American Institute of Aeronautics and Astronautics 
and the Confederation of European Aerospace Societies 
Breckenridge, Colorado, June 17-19, 2002 


National Aeronautics and 
Space Administration 


Glenn Research Center 


July 2002 



Acknowledgments 


The authors would like to thank Dr. Stewart Leib, NASA Glenn Research Center, for helpful comments in writing 
the manuscript. In addition, the authors wish to thank Dr. Srini Bhat, Boeing, for providing the jet noise data. 


The Aerospace Propulsion and Power Program at 
NASA Glenn Research Center sponsored this work. 


Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 


Available electronically at http: / / glti~s.grc.nasa.gov /CL.TRS 




A PARAMETRIC STUDY OF FINE-SCALE TURBULENCE MIXING NOISE 


Abbas Khavaran 
QSS Group, Inc. 
Cleveland, Ohio 


James Bridges 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 


Abstract 

The present paper is a study of aerodynamic noise spectra 
from model functions that describe the source. The study 
is motivated by the need to improve the spectral shape of 
the MGBK jet noise prediction methodology at high 
frequency. The predicted spectral shape usually appears 
less broadband than measurements and faster decaying at 
high frequency. Theoretical representation of the source is 
based on Lilley’s equation. Numerical simulations of high- 
speed subsonic jets as well as some recent turbulence 
measurements reveal a number of interesting statistical 
properties of turbulence correlation functions that may 
have a bearing on radiated noise. These studies indicate 
that an exponential spatial function may be a more 
appropriate representation of a two-point correlation 
compared to its Gaussian counterpart. The effect of source 
non-compactness on spectral shape is discussed. It is 
shown that source non-compactness could well be the 
differentiating factor between the Gaussian and exponential 
model functions. In particular, the fall-off of the noise 
spectra at high frequency is studied and it is shown that a 
non-compact source with an exponential model function 
results in a broader spectrum and better agreement with 
data. A recent source model proposed by Tam et al. that 
represents the source as a covariance of the convective 
derivative of fine-scale turbulence kinetic energy is also 
examined. 

1. Introduction 

A physics-based methodology utilizing the averaged 
equations of motion is used to assess jet noise spectra as a 
function of the source model function. The focus here is on 
small-scale turbulence noise that dominates the spectra at 
subsonic Mach numbers. It is generally accepted that 
sound generation in jets is a by-product of the unsteady 
features of the flow. Any flow manipulation intended to 
reshape the spectra and subsequent perceived noise level 
ought to affect the unsteady characteristics of the flow. 

Mixing enhancement devices such as chevrons and tabs are 
known to alter turbulence statistics 1 and change the time- 
and length-scales of noise generating eddies. These flow 
modifications directly impact the sound field by, for 
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example, reducing the low-frequency noise at the cost of 
adding to the high-frequency content. 

As such, it is not unreasonable to argue that modeling of the 
unsteady behavior of the flow as supported by a number of 
isolated test configurations may not lead to a reliable 
prediction tool. In an ideal situation, one might attempt to 
solve the full compressible Navier Stokes equations without 
resorting to any modeling closures. In fact, it is simply a 
matter of time before Direct Numerical Simulations (DNS) 
should solve the far-field jet acoustics either directly or by 
some extension of the near-field solution 2 . Clearly, the 
computational demands of a typical high Reynolds number 
jet leave little room, at least in the near future, for DNS as a 
design code. Nevertheless, it could be an extremely useful 
tool in understanding the unsteady features of the jet in order 
to improve the source modeling. 

On the other hand, the computational requirements are 
greatly reduced if the governing equations are spatially 
filtered, as is done in the Large Eddy Simulation (LES) [e.g., 
Ref. 3, 4, 5], and the effect of subgrid scales (SGS) is 
modeled. This approach is successful in capturing the 
distinct directivity of supersonic jets that results from large 
scales of motion (instability waves) and dominates the 
general noise picture near the down-stream axis. But it might 
also suffer from a neglect of high-frequency noise resulting 
from subgrid scales. A recent study by Seror et al 6 calculates 
the acoustic pressure from the filtered Lighthill’s stress tensor 
as well as the full tensor and concludes that the SGS part of 
the tensor needs to be taken into account in order to recover 
reliable high-frequency results. One must also be careful 
about the selection of the subgrid scale eddy viscosity model 
e.g., constant vs. dynamic Smagorinsky coefficient, as well as 
the inflow perturbations used in the simulation. LES 
predictions of Morris et al. 3 appear to benefit from the 
dynamic model and narrow down some of the usual over- 
predictions in turbulence and noise. 

Physics-based prediction methods such as MGBK 7 or Tam 
and Auriault’s 8 fine-scale model heavily rely on model 
functions that express the statistical properties of noise 
sources. These predictions use the averaged equations of 
motion; hence the unsteady features of the flow are entirely 
described by two-point, space-time correlation models. Any 
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shortcoming in the predictions should be directly linked to 
the model. Other phenomena such as refraction and 
convection impact the directivity of jet noise. 

Woodruff el al. 9 study the isotropic source model in the 
original MGBK noise prediction methodology 10 and 
examine alternative representations for the turbulence 
spectra. They propose an energy spectrum function for the 
two-point velocity correlation to satisfy the Kolmogoroff 1 1 
spectrum law in the inertial sub-range. In particular, they 
examine a wave-number-dependent Gaussian function for 
the temporal part of the correlation, as opposed to the usual 
separable space and time functions. Here the characteristic 
frequency is scaled based on the spatial wave-number k and 
the turbulence dissipation rate £ as t ~ k J 'e' L \ The 
predicted spectra do not appear to offer noticeable 
improvement over the conventional MGBK methodology 
that uses a separable correlation function and calculates a 
characteristic frequency from the turbulence kinetic energy 
k as t” 1 ~ k / e . Nonetheless, near 90° angle, they report 
slight improvement in their spectral shape. 

Tam et al & model their source as a two-point correlation of 
the convective derivative of kinetic energy of small-scale 
turbulence. They use RANS to calculate the time- and 
length-scales of the noise sources as is done in the MGBK, 
and predict noise spectra in good agreement with data at 
mid angles. 

In all, the noise from small scales of motion, which are 
usually broadband in nature and cover a range of observer 
angles, remains a significant part of the jet noise spectra. 
Indications are that at high subsonic Mach numbers, and 
heated jets in particular, instability-associated noise may 
dominate the low end of the spectra at shallow angles. On 
the other hand, one might argue that the mean-flow effects 
could also play a role, by diverting the high-frequency 
noise of small-scale turbulence away from the axis and 
creating a region near the zone of silence that is dominated 
by low-frequency noise. 

In this work, alternative representations of the source in 
modeling the quadrupole terms of Lilley's equation are 
examined. We concentrate on 90° emission angle where 
shear-noise is not a factor and the spectral shape is 
predominantly defined by the self-noise component. 

The paper begins with some preliminary formulation of the 
governing equations. Section 2.2 derives expressions for 
source and non-compactness factor as a function of the 
proposed models. It is shown that the fall-off of the high- 
frequency noise becomes less steep when an exponential 
spatial function with an appropriate temporal function is 
selected and non-compactness effects are included. Section 
3 compares Tam and Auriault's model 8 with the MGBK 
model as proposed in Ref [7]. It concludes that the spectral 


shape at 90° should be identical if consistent assumptions are 
used. Some concluding remarks on future directions for a 
physics-based modeling approach is given in the summary. 

2, Sound Spectral Density 

Application of Lilley’s equation to the problem of jet noise 
and the significance of various source terms has been the 
subject of numerous discussions in aeroacoustics. In a recent 
article, Goldstein 12 gives an exact form of the equation with 
the dependent variable defined such that the source is of 
quadrupole/dipole nature. The quadrupole source is second- 
order in velocity fluctuations and is the sum of the commonly 
known self- and shear noise terms. The dipole term is 
produced by the fluctuating sound speed due to temperature 
fluctuations. 

The far-field spectral density due to sources of Lilley's 
equation may be expressed as integration over the source volume 

y 

P ( X I o» = J j G\x, y - 1 / 2, co)G{x, y + \ / 2, (o)Q n (y, \ a) 

i ( 

xd^dy. (1.1) 

G is an appropriate Green’s function, * denotes a complex 
conjugate, and Q n is the source spectral density which is 
formed from a Fourier transform of a two-point space-time 

correlation between source points y, =y-%/2 and 
y, = y + 1 / 2 separated by time T 

+oo 

Q u (y,lw) = jR(y,lt)e im dt . (1.2) 

If we assume that the variation of the magnitude of the 
Green's function with respect to within the source region 
where Qn is non-zero, is negligible compared to that of its 
phase then 

G\x,y / 2 ,ot)G(x,y + 2,(0) ~ |G(x,y,ffl)| 2 e~ ,k *, (1.3) 

where k is a wave number of magnitude (c o la _) and direction 
(x-y). Here co is the frequency at tire observer location and a _ 
denotes tire ambient sound speed. Equation ( 1 . 1 ) is no w written as 

p 2 (x\w) = j\G(x,y,w)\ J e'^dtj R(y,^,t)e~‘ ki d^dy. (1.4) 

j i 

It is argued that changes in retarded time across a correlation volume 
element are more likely to be small if the correlation is written in a 
frame x ' moving with convection velocity U c (i.e. x' = x - i U c t ) 

of tire turbulent eddies. Experimentally, the correlation R(y, §, T)in 
a jet flow describes a fluctuating pattern in a moving frame and is 
expressed as 
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R(y,l,r) = RJy,l m , r), l m =l~iU c r. 


(1.5) 


C = jg(r,d)dr. 


(2.4) 


Source frequency Q is related to the observer frequency 
through the usual Doppler effect Q = o(l — M c cos 6). 
Making a transformation to the moving frame, we find 
-k. \ + cm = -k. + Qt, therefore 

R(y,^,z)e~' k4 d^=je aT dzj R m {y,^ m ,z)e~‘ k ^’d^ in 
« — L 

(1.6) 

Now the noise spectral density with respect to the moving frame 
is given as 




p 2 (x I (o) = J | G(x, y , Q) |" J e iat dz J R m (y,% m , r)e 

y ” S. 

(1.7) 


2,1 Self Noise Spectra 

For brevity, subscript m is suppressed in this section, and § is 
used as the separation vector with respect to the moving frame. 

To assess the far-field mean-square pressure due to the self-noise 
term only, we write Lilley's equation in a coordinate x' moving 
with convection velocity U c 


Subscript o refers to source location and the shielding function 
g( r , 6 ) is defined in Appendix A. 

The above expression for S is applicable outside the zone of silence 
of a source only, where the shielding function g“ (r, 6) is positive at 
all radial positions; hence there is no shielding. The acoustic 
pressure due to the above source and Green's function becomes 

d\u.u.) 

P, ef (x ,0=11 5 , (x',l;y,l 1 )p ' J dt,dy, (2.5) 

L y Jy.dy : 

S{x ,t\y,t j ) is the inverse Fourier transform of Six' , y,&>) 

1 +oo 

S(x ,t;y,t . ) = — fs(jr t y,£2)e"“ l '" , ’dn (2.6) 

2;r i 

Upon transferring the derivatives from the source to the Green’s 
function in (2.5), and making the approximation that the variation of 
the magnitude of the Green's function Six ' , y, co) with respect to 

separation vector £ within the source region is negligible compared 
to that of its phase, we find 


L(p sel /,V,x l )= D 


( d 2 (u i u) 2 ' 

P 


V 


dx.dx. , 

1 J / 


( 2 . 1 ) 


Here L is Lilley operator, V=U - U c , and die density p has been 
moved to the right of operator D = 3 / dt + Vd / assuming 
that flow is locally parallel and that density fluctuations are small 
so p is the mean density. The Green's function to the above 
equation for a source of type D{e‘ n, 8{x - x o )}is 

USe^-V,!^ = D{e" n, S(x 1 )8{(p-(pjS(r-r)/ r}. (2.2) 
In the high frequency limit, S is given as 

1 ala 


Pi,/ (* I ®) = J I (x' , y, D.)S, kl (x' , y , H)l I ijU (y , Q)dy . 
y 

(2.7) 

Subscripts on S refer to derivatives with respect to source 
coordinate y, and the phase factor e~‘ k 4 is now included with 
the source correlation I iJkl . In addition, we have neglected 

the mean density gradients so that density p is included 
within source correlation function 


S(x ,y,Q) = ■ 


4 nR (1 - M c cos 0)(1 - M o cos d) 
Q 




= P~j J (UiUyX^u,) e~' k4 e ni drdl . ( 2 . 8 ) 

i 


The volume integration in (2.7) usually includes the most energetic 
parts of the jet. For axisymmetric jets, the directivity factor may be 
averaged azimuthally with respect to source and observer 
circumferential angles to obtain a ring-source directivity factor, 
a jjkl . Subsequently, jet volume integration will be limited to radial 

and axial coordinates 


X e iau *- exp{/ — [ j ; (g - gJdr-Ccos((p -yj]} (2.3) 

a 


*ijkl 


-y + 71+71 

= j j I ^ , i j k l \ d(pd(p a , 


(2.9) 


where 


p 2 (x\a>) = J J ^jkikjkeiS^ QY'lnrdrdyJ. (2.10) 


— >0 
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Equation (2.10) is written in an expanded form for a unit ring- 
source volume at point y 


X ( Sel f N ° ise > = 7 mi a n 


7 2222 ^ 2222 


7 3333^3333 _ 7 122 'L 122 


" 77 I2I2^I212^ “ 7 | 1 3:7 1 I3:7^ 7 |3i:7hi:7 “ 7 2233 ^2233^" “^ 7 2323 ^ 2323 ’ 


( 2 . 11 ) 

Directivity factors a. jkl , and shielding coefficients p are 
defined in Appendix A. In axisymmetric jets. 


amounts to setting the wave number equal to zero. Thus the four- 
dimensional transform (2.8) simplifies. 


I ijkl {y,il) — P 2 J J («,« )(u k u , ) e^dTdl 

i — 

(compact-eddy). (2.15) 

Using an appropriate model in (2.13) with L = 0 one finds 


^2222 ^ 3333 ’ ^1 212 ^1313 7 122 7 133 ’ 

^2233 ““ ^2323 

At 90° (2.11) simplifies 


a -Am(y.n) 


^=p 2 («;) 2 f 3 G(Q), 

— p 2 ( Ml 2 ) 2 f 3 G(D), 

5tt 


£ 2 

/(£) = exp(-7r— ) 

f 

/(I) = exp(-TT-) 

£ 


X (SW/ Atofre) = 2/ 2222 a 2222 + 2(/ 2233 + 2 / 2323 )a 2233 , 

0 = 90° . (2.12) 

At this point, a physics-based modeling approach is employed to 
obtain closed-form expressions for the correlation coefficients 
appearing in (2.11) and (2.12). For convenience, we assume that 
the joint probability distribution of velocities u and u’ at points 
y and y ' (separated in space and time) is normal and write the 
fourth-order cross-correlation function as a superposition of 
second-order correlations B . In addition, the second-order 

correlation is assumed separable, i.e„ u i u j = R ij (^)g(t) . 

Here /? (£) and g(r) denote the spatial and temporal parts of 

the correlation, respectively. From (2.8), the axial correlation 
coefficient becomes 


/ im (y,Q) = 2p 2 G(Q)J<(^-^ dl 

i 


+oo 

G(Q) = j g 2 (T)e ar dT. 


(2.13) 


In homogeneous isotropic turbulence, the two-point 
correlation function has the form 1 1 


Rijti) = u;[(f+^f')5 iJ /& ■ (2.14) 


< g2 g2 

c. / im (y,«) = — r p\u l ) 2 f./ 2 G(Q), r(|) = exp{-;r(% + %} 
2V2 £ £ 


d ■ / llu (y,£2) = — -P(KY ( A G ( Q )> r (<^) = exp{-^ 


5 n 



with ^ 2 a =^\+^\. Models (a) and (b) assume isotropic 
turbulence, whereas (c) and (d) use an axisymmetric 
turbulence constructed from a set of kinematically compatible 
scalar functions [Ref .7], Length-scales f 3 and £ 2 are 

2 3/2 2 3/2 

proportional to (u l )' 7 e and (m 2 ) / e , respectively. 

Consider the isotropic model, cases (a) and (b). After 
substituting the remaining correlation coefficients into (2.12) 
(see Appendix B) we find 

^(Self Noise) = 2 (a 2222 + a 2233 )/ ml . 9 = 90° (2.17) 


The shielding coefficients /D that multiply directivity factor 
a ... result in noise attenuation within the zone of silence. At 

ijkl 

90° emission angle, Q = to , and the shielding function is 
g'(r) = (a m / a) 4 , where a is the sound speed at the source. 
Here g 2 (r) has no zero crossing, therefore there is no turning 
point and /1 ; = 1 . 


In the following discussion we examine Gaussian as well 
as exponential spatial functions /(£). 


Y (Self Noise) =ISf k\—) 4 I uu , 

a 


6= 90° (2.18) 


2,2 Compact Eddv Approximation 

In a compact eddy approximation, the assumption is made that 
the eddy length-scale I is much shorter than the wavelength of the 
acoustic disturbances, i.e. all is small compared to unity. As 

such the factor e~' k 4 is set equal to unity, which practically 


Away from 90° equation (2.11) should be used directly. In 
addition, the shear noise contribution needs to be included as 
described in Ref [7]. Model (c) shows that the directivity of 
jet noise with respect to its level at 90° is a function of 
anisotropy of turbulence 7 ' 
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Let t be the inverse of the characteristic source frequency, 
which is proportional to turbulence kinetic energy and its 
dissipation rate as t = Q = a i e / k . Eddy length-scale is 
obtained from l ~ % u x . At this point noise spectrum at 90° 
may be calculated using models (a) or (b) in (2.18) and 
with an appropriately selected proportionality factor for i . 
It is clear that the spectral shape will be the same for both 
models. For instance, if the proportionality factor selected 
for length-scale f. in model (b) was larger than that of 
model (a), say by factor (5 jt‘ /8-\/2) 1/3 , then identical 
spectra are obtained. 


Upon examining the three-dimensional energy spectrum for 

a two-point correlation R tJ {£) using models (a) and (b), it 

is found that both scale as fourth power of spatial wave 
number when wavelength is large. However, in the inertial 
sub-range, their decay rate is substantially different. The 
energy spectrum function, normalized with respect to 

(1.5 u[), is given as 


E(k) = 


f (, kf. f k 2 f 2 

pexp(- ), 

12 n {Arif An 


E(k) = 


16 f {Ml n) 


3ri 


k 2 f 2 


(i + — r) 


/(|) = exp(-jr— ) 

f 

£ 

/(|) = exp(-jr-). 

e. 


(2.19) 


Figure 1 shows that the energy spectrum for the 
exponential function (model b) decays as k at large 
spatial wave number. This is reasonably close to the 
proposed Kolmogoroff's A’ ~ 5/3 scaling law 11 . Model (a), 
however, presents a much faster decay rate. These 
comparisons suggest that when non-compactness effects 
are accounted for, the two models, if properly used, should 
produce different spectral decay at high frequency (this will 
be shown in section 2.4). 


Reference [11] argues that the exponential function may 
not be strictly correct on the grounds that (1) it is not 

parabolic at its vertex, (2) the lateral correlation with 

£ in direction of £ 2 remains positive for all £p, whereas 
the correct curve must become negative for large £ 2 . 

Recent measurements of Bridges et al . 14 appear to suggest 
that a two-point correlation R tj {£) constructed from an 

exponential spatial function according to (2.14) provides a 
better fit to data relative to the Gaussian function (Fig. 2). 


2,3 Source Non-Compactness 

Here we explore the effect of source non-compactness on 
noise spectra using Gaussian and exponential models (a) and 

(b). To carry out the integration with respect to £ in (2.8) or 
(2.13), a convenient coordinate transformation has one of the 
axes £ { aligned with vector k . 

In a spherical coordinate £ = £{cosa, sin a cos^t, 
sin a simp), £ x is aligned with wave number k such that 
k.£=k £ cos a . Equivalently, one might select a cylindrical 
coordinate system £-(£ l , rcoscp, r sirup) with £ x in the 

direction of k such that k.£=k£ x . In any event, when 
turbulence is isotropic, the final result should be independent 
of the direction of wave number k . Now 
with/(£) = exp {-n£ 2 If 2 ) we find 


A.iiCu^) 


2VI 


p 2 {u 2 f f 3 G{Q)N{kf ) , N{M) = exp(- 


k t 

) 

8 n 


( 2 . 20 ) 


where N{kf) denotes the non-compactness factor. Using an 
exponential function / {£) = exp(-^ / (.) we find 

I uu {y,Q) = —p 2 {u 2 ffG{kl)N{kl), 

Sri 


k(. 2 

n 5 ! k(. kf. 5( ^ } +12 

N{k(.) = 20( — ) [3tan ( — )-2 

kf 2n n 


kf 2 

(— ) +4 
n 


( 2 . 21 ) 


Applying the law of limits repeatedly, the last expression 
iorNfkf.) becomes 1 as kf approaches zero. Figure 3 shows 
that the non-compactness factors equal 1.0 for 0< kf< 2 and 
decay rapidly for kf>l . The Gaussian function appears to 
produce a faster-decaying non-compactness factor. The 
above non-compactness factors repeat for other correlation 
coefficients as noted in Appendix B. 

Since f ~ T o «,, we find kf ~ (c QT o )(u l / a m ). Ratio u x / a„ 
is usually smaller than 1.0 (of the order of 0.2 for the more 
energetic parts of the jet). In the compact eddy 
approximation kf was assumed small henceN{kf.) was set 
equal to 1.0 for the entire range of the wave number. As 
frequency CO becomes very large, (OT o may be large enough 
to reduce A below 1.0 (Fig. 3). 

The implications are that the effect of N{kf) on spectra, if 
any, should be a slower decay rate at high frequency for the 
exponential function relative to the Gaussian. It is also 
evident that both models produce a faster decay at high 
frequency when the sources become non-compact. However, 
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as we shall see shortly, these effects are visible only in the 
context of the spectral shape function G(Ll)N(kl). If 
G(Q) has already decayed far enough before N(ki) takes 
effect, then both models produce identical spectra. 

Figure 4 shows the MGBK prediction using Gaussian- 
isotropic source model (a). The effect of source non- 
compactness on predicted noise spectra for a Mach 0.5 cold 
jet is of the order of 0.12 dB at the high end of the spectra. 
Here the temporal part of the correlation was selected as 7 

g(T) = exp{— ^(<7 / 2) 2 +(r / t o ) 2 } . (2.22) 


Fig. 7 is the MGBK predicted spectrum for Mach 0.5 cold jet 
with <7 = 0, including the non-compactness. Model (b) 
predicts a broader spectrum and noticeable high-frequency 
improvement relative to model (a). Here, the location of the 
peak frequency was adjusted slightly by selecting 
proportionality constant a l , =a l £/K) as 0.225 and 
0.170 for models (a) and (b) respectively. 

3. Tam and Auriault's Model 

In reference [8], Tamet al. compute the fine-scale turbulence 
noise from an equation similar to (1.4). Here we compare the 
Green's function as well as the source cross correlation 
functions between the MGBK model and Tam's approach. 


Constant <7 = 0.8 as was originally proposed in Ref [7], 

One might expect a similar effect at other angles, 3.1 Green s Function — 

asN(kf.) is a common factor throughout Eq. (2.11). The Green s function to the linearized Euler equations for a 

locally parallel flow is the solution to 


Next we explore the MGBK spectral shape function. 


2,4 Spectral Shape Function 

Apart from factor k 4 , which appears due to transfer of 
derivatives from source to the Green’s function (not 
included in the following discussion), the spectral shape 
function, denoted as F^onJ, is simply the product 
G{kl)N(ki) . Function G(Q) is obtained from (2.13) and 


( 2 . 22 ) 


= err. 


K t [o-yjl + (Qt o / 2) 2 ] 
Ji + (QtJ2) 2 


N{k(.) 


(2.23) 


The normalized spectral function becomes 


L(Ge~ ia '\U,x,) = e~ im, 8{x-x s ) , (3.1) 

where CO denotes source frequency with respect to stationary 
frame X, x s is the source location, andL is Lilley's operator. 
Tam and Auriault 15 recast the problem into an adjoint 
operator for the adjoint Green's function G a , which is related 
to the Green's function G of the original problem by a simple 
switch of the source and observer locations X,. and x„ . The 
final result for an axisymmetric mean flow (now multiplied 
by 2k(C to compare with the high-frequency solution) is 
given as 




-y/l + (Qr, / 2) 2 K^a) 


N(kf), 


Q = o(l- Af. cos0) . 


G(x 0 ,x s ,co) = G a (x s ,x 0 ,co) 


-ik(xcosQ-R) 


(2.24) G a (x s ,x 0 ,CO) 


Anka R 




cos m<p ■ 


r<R. 

s o 

(3.2) 


As <7 becomes very small one finds 


C.im F x 

(7 — >0 


' :A W) 

i+lQzjiy 


Function f m (r) is obtained by solving an ordinary 
differential equation and matching the above inner solution 
(2.25) with the outer solution at the jet boundary^ . 


Figure 5 shows the spectral shape function F x (cor o ) for 
model (a), with t. /(UT o ) = 0.20, which applies to the 
more energetic parts of the flow. Fig 5a shows a negligible 
role for the non-compactness factor at <7 = 0.80 . The 
effect of N{kk) becomes evident as shown in Figures 5b 
and 5c. It should be noted that when <7 — > 0 the compact 
source model becomes increasingly inadequate (resulting in 
unusually high level of noise at high frequency). 

Comparison of the spectral shape functions of models (a) 
and (b) shows a slower high- frequency decay for the 
exponential model (Figs. 6a and 6b). Shown in 


-ik(xcosQ-R) of 

G a (x,x o ,( 0 ) = 6 X (Z) +AX 11 <*)] 

4?z -ka„R m=0 

xcos mep, r>R o (3.3) 

where % ~ kR o sin 6 . 

The corresponding high-frequency solution 16 for a non- 
convecting (£2 = co) monopole type source takes one of the 
following forms depending on the location of the source r o 
relative to the zero crossing point r a of the shielding function, 
i.e. g(r a ) = 0 (see Appendix A for the definition of g). 
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G(x,y,co) = 


a la 




AnaJtR (1- M o cos 6) \rh o J 
X exp( ik[\; (g - )dr - gj a ] - k[£ a - C 0 cos «p - <p o )]), 


r <r 


(3.4) 


with h = -g , 
or 


£ = J h(r,6)dr 


G(x,y,co) = - 


i 


ala 


4 najiR (1 - M cos 6)~ 


/ \ 1/2 

' £ A 

o 




xexp(ik[\~(g-gJdr-C 0 cos((p-(pJ]), r > r 
with £ = j^g(r,6)dr. 


(3.5) 


section we use £ as separation vector in a fixed reference 
frame. 

The corresponding axial correlation coefficient used in the 
MGBK describes the cross correlation of the Reynolds stress 
components rather than their convective derivatives. 
Following the usual MGBK methodology, the fourth-order 
correlation / 11 u (£,t) = ( i (pu 1 u 1 )(pu l u l )j is expressed as a 
sum of second-order tensors. The element of I nn 

contributing to the noise field is simply 27? [ 2 [ (£) g 2 (T) 
which, upon using model (a) in (2.14), and making a 
transition to a fixed reference frame, becomes 

/ u„ (l T) = 2 P 2 (u 2 y (1 - J £; ) 2 exp{ - 2 ^ J [(£, - Ur) 2 

+ £ 2 3 ]}g 2 (T) (3.7) 


Figure 8 shows comparisons of the above Green’s 
functions for a stationary ring source at 7 diameters from 
the jet exit ( x/D = 7), and with radius r o / D indicated as 
a parameter. Computations were carried out for a range of 
Strouhal numbers (St = fDI LF) for a Mach 0.9 cold jet. 

The agreement is generally good at high frequency. As the 
Strouhal number is lowered, the high- frequency 
approximation appears to deteriorate. Discrepancies 
become increasingly visible near the boundary of zone of 
silence as seen in figure 8c, with the adjoint Green’s 
function predicting a larger zone of silence. For our 
purpose, we intend to compare the MGBK spectrum with 
Tam’s solution at 0 = 90° and assess the high-frequency 
behavior. Spectral peak for both Mach 0.5 and 0.9 cold jets 
is near St = 0.9, therefore the high-frequency solution is 
found suitable. 

It is noted that the Green’s function for a monopole type 
source scales as A'* 1 with respect to the wave number and 
as 1 / (1 — M o cos 6 y with respect to polar angle. 


and temporal function is given by (2.22). 


Factor (l-^£ 2 3 if) in (3.7) reduces / u11 (£,t) to zero as 
the normalized lateral distance (£ 23 If.) approaches 1/ yfn . 
Beyond this point the correlation is practically zero. 

In order to compare correlation functions (3. 6) and (3.7) on 
an equal basis, suppose we relate time- and length scales and 
define the following dimensionless parameters 

2 ji — 

f. = (— — ) 05 f , T=2 t, !.=!•//, T = T/T . 

fnl ’ ° ’ ‘ ‘ ’ 

(3.8) 

Now the lateral correlations, with zero time-delay, and 
normalized in magnitude become 


/ Dqfxjy Dg,(j + £ 2 3,g \ 
\ Df Df / 


exp{-7« 2 £ 2 2 3 } 


(3.9) 


3.2 Source Model 

Tam and Auriault 8 propose a two-point, fourth-order, axial 
velocity correlation in a fixed reference frame 

Dqfx^ty Dq s (x 2 ,t 2 )\ q s r -||J 

Dt Dt , / ex U t 

1 2 Is s 

f.n 2 , , 

- — [(^-Uxy+O) (3-6) 

with £ = X, - X 2 , X = t l —t 2 , £ 23 as defined earlier, and 

U the mean velocity at the source location, which can 
be replaced with the source convection velocity. In this 


and 

fn2— — 

W5 * >0) = (1 ~ O' exp{-^ 2 £ 2 3 } . (3.10) 

2 

Figure 9a shows that (3.9) and (3.10) decay somewhat 
differently. Tam and Auriault’s correlation (3.9) does not 
have a zero intersect and decays at a slower rate. Limited data 
available on fourth-order correlation measurements 17 seem in 
better agreement with the MGBK model. However, this 
difference may practically be insignificant in noise 
prediction. An interested reader may find more information 
on the second-order lateral correlation in a book by 
Townsend 18 . 
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Aside from <§ 23 , which is now set equal to 0.0, the two 

models appear similar. The normalized axial cross 
correlation functions are: 


j Dq z (x,t j) Dq s (x+ ^,t 2 ) 

\ D h Dt 2 

-In 2^ 


= exp{-a|^| 
r / or) 2 }. 


(3.11) 


Am (<?! , T) = expj-V<T : +T 2 -/j»2(| 1 -T/a) 2 } 
with a =(f. s / Urj. 


(3.12) 


Following Eq. 33 [Ref. 8 ], Tam and Auriault's spectral 
function F, , is written for a unit volume of turbulence at V 

— I £ I 

F 2 (ar ! ) = f f p a (y,x,-co)p a (y + l,x,+(D) rf|exp{- — 

t ; 

-—[(§! -t/T) 2 +^] + ifl)T}rfrd|. (3.15) 

Here p a is the direct Green’s function, which is obtained 
from switching source and observer locations in the adjoint 
problem. As was done earlier, the product of the Green's 
functions is approximated as the magnitude at the center of 
the correlation volume multiplied by a proper phase 


Using Tam’s constants (£ s = c , K 312 / £, T s = c t k / £, 
c f = 0.256, c z = 0.233) one finds a = (c f / c t )(k 05 / U). 
Let’s set U equal to the convection velocity 0.65 U } . In the 
more energetic parts of the flow (mixing layer) we chose 
( K°' s l Uj) = 0.12, and find a = 0.20. 

Figure 9b shows Tam’s correlation coefficient (3.11) and 
the MGBK model (3.12), with O = 0. The effect of small 
parameter O is shown in figure 9c. Aside from the slight 
difference described above in comparing the lateral 
correlations, the two models exhibit similar features. 


p a (y,x ,-(o)p a (y + l,x,+co) » \p a {y,x,-cof e " 4 . (3.16) 

Phase factor k . £ should not carry a preferred direction since 

rays are emitted at all angles and separation vector £ may 
also take any direction. It follows 

F 2 (m,) = \p a (y,x-(o )\ 2 J Jexp{ - --^[(l, - Ur ) 2 

£ T U T S f‘s 

+ t; 2 23 ] + ik.%+iar}dr , 

or upon integrating over T 


Next, the spectral shape functions are compared. 

3.3. Spectral Function 

The MGBK spectral function was written earlier in a 
moving frame as a Fourier transform of the temporal 

function g"(T) , multiplied by the non-compactness factor. 
In a fixed frame variable T appears in (|, - Ur ), hence 
the spatial function needs to be include in the integration. 
However, with a simple transformation = E,-Ur (see 
1.6), the spectral shape function (2.24) is recovered. For 
comparison with Tam’s spectra we now let T g = 2t s and 

relate t to f. s (see 3.8). 

[aJl + (QrJ] 

v N(k(), 

L + (£2t s ) 2 K 1 ((x) (3.13) 

Q = co(l-M c cos 6). 

N(M) for model (a) is given in (2.20). As a become 
small one finds 
1 

tim F = 7 N(k() . (3. 14) 

°->° 1 + (£2 t s ) 



F 2 (cot ) = \p a (y, x-cof expf — (— ) 2 ) J exp{ - 

1 1 /?., "I 'ITT J TT~ 


in 2 2 U 


Ur 


(m2 ©I _ 5 

%*~ l ~ + lk -^ d S- 

’ s 

As before, make a coordinate transformation with respect to 
dummy variable <fj such that aligns with k (i.e. 

k.E, = k^). The component of U in direction of k becomes 
(A = U cos 0 . 


F 2 (wrj=\p a (y,x,-w)\ 2 2f 3 s expf- — (— ) 2 If J^ 2 rf | 3 

\ (m2 2U J 


: J exp{- 


(E, , co 

-fn 2 | 23 }cos[( k)4jd^ . 

o U e r s _ U e 

Next transition to polar coordinates ! = (!,, r cos 9 , r sirup) 

2 nf Ur K. 


F 2 («n s ) = \p a (y,x-G))\ 2 


(m2 1 + (wr -U 0 rk) 


xexp 


1 co( 

— (— ) , 
\ (m2 2U J 


(3.17) 
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Expression (3.17) is now normalized to obtain the spectral 
shape function. We also take out the magnitude of the 
Green’s from the spectral function (as was done in 3.13). 
For a compact eddy ( k = 0) 


F 2 {on,) = 


1 

1 + (<» T ,) 2 



1 (Of 2 

(— ) 

In 2 2 U 


A 


(compact eddy) (3.18) 


obtained by replacing the k 4 wave-number factor by 1C (as 
in Tam and Auriault approach) are also shown. In doing so, 
some minor adjustments had to be made in the calibration 
constant related to the source characteristic frequency to 
persevere the location of the peak spectra. The noise 

spectrum naturally becomes broader with the ZT factor 
(compare 1 and 2), and improvements are noticed at both 
ends. 


and 

F 2 (cotJ = 


1 


U 


■exp 


^ 1 cot 2 ^ 

(— ) 


l + (<MT s ) (1 COS 6 ) 

a 


V. f.n 2 2 U J 
(non-compact eddy). (3.19) 


This result is the same as that reported by Tam and 
Auriault, but was derived without resorting to the following 
approximation, which was suggested in Ref. [8], 

p a (y,x-(D)p a (y + |,x,+<») = \p a (y,x,-co)\ 2 exp {/*£, cos 0} 


At this point let’s compare (3.19) with the MGBK shape 
function (model a) at 90° emission angle. Figure 10 shows 
comparisons with l s I Ut s = 0.20 as was selected earlier. 
The MGBK spectral shape (with a = 0.20 ) agrees quite 
favorably with Tam's results at U / a„ =0.325 and 0.65 
as shown in figure 10a and 10b, corresponding to jet exit 
velocity of U ; I = 0.50 and U } la,. = 1.0 

respectively. 

It is noted that in comparing the above spectral shape 
functions, we have deliberately dropped the wave number 

power k 4 from MGBK model (see 2.18), and k 2 from 

Tam's model. In the MGBK approach, factor k 4 appears 
when two spatial derivatives and one convective derivative 
are transferred from source to the Green's function. Tam 
and Auriault, on the other hand, maintain that convective 
derivative is included in their source modeling. As a result 
the Green's function p a remains of the order of k (see 
equation (28) of Ref [15]), which incidentally indicates that 
source q s should be of quadrupole type. With the 
convective derivative now hidden in the source, the power 

spectral function becomes proportional to k \ 

Morris and Farassat 19 described this in more detail in a 
recent paper and suggested that a consistent approach 

should result in ~ k 4 for both MGBK and Tam's analysis. 

Figure 11 shows the MGBK predictions for Mach 0.5 cold 
jet using model (a) and a spectral shape function that 
matches that of Tam and Auriault's (curve 1). Predictions 


Additional high-frequency improvement could be gained by 
simply removing the atmospheric attenuation built in the 
MGBK code from predictions and implying that it is built 
into the source model (as was done in spectral predictions of 
Ref. 8 at R/D = 100). However, predicted noise spectrum 
should account for atmospheric attenuation. This amounts to 
attenuating the high-frequency noise depending on the 
observer distance and atmospheric conditions (i.e., relative 
humidity and ambient temperature). Noise measurements 
usually reflect the atmospheric attenuation. Figure 11 (curve 
3) shows the significance of atmospheric loss on predicted 
spectra. Although excellent agreement with data is thereby 
obtained, two very questionable steps were taken to predict a 
better spectrum. 

As was shown earlier (Fig. 7), good agreement with data 
could be achieved by selecting an exponential spatial function 
(model b) in place of the Gaussian function. 

4. Concluding Remarks 

In the preceding discussions, we examined alternative model 
representations for the two-point space-time correlation 
appearing in physics-based jet noise prediction 
methodologies. It was argued that a proper representation of 
the source, consistent with the observations and accepted 
fundamentals related to turbulence statistics, should improve 
prediction of the flow-generated noise in the framework of 
Lilley's equation. The discussions centered on noise from 
small-scales of motion and at 90° observer-angle. The main 
result of the study was described in Fig. 7. It was shown that 
an exponential spatial function, with source non-compactness 
included, predicts a broader spectrum relative to a Gaussian 
function. The effect of turbulence anisotropy may readily be 
accounted for by selecting model (d) of section 2.2 as a non- 
compact source. 

Away from 90°, mean-flow refraction effects as well as 
convective amplification due to source motion become 
crucial in capturing the peak directivity that occurs near the 
down-stream axis. A high-frequency approximation offers an 
analytical solution to the Green’s function, but comparisons 
of section 3.1 indicate that it might not be an appropriate 
approximation at small Strouhal numbers. 

A numerically computed Green’s function 15 provides extra 
flexibility at low frequency, at the cost of added numerical 
intensity. 


NAS A/TM— 2002-21 1696 


9 



Appendix A 

Directivity factor for various quadrupole source 
components (in the absence of mean density gradient) is 

cos 4 6 Q , 2 n 

v=— — — (— ) IS1 2 p n , 

(1-M c cos0) a _ 
g 2 (r) cos 2 6 Q 4 , 

fl 1122 = ° -(— ) 4 \s \ 2 p u , 

2(1 - M c cos6) a. 
fl 2 2 22 = (3/8)g 4 (r )(— -) 4 ISf P 22 , 

fl 22 33 =(l/8)^ 4 (r)(— ) 4 I5| 2 j8 23 . 

The shielding function is 

(1-M cos 6) 2 (a ^ / a) 2 -cos 2 6 


8 (r,6 ) = 


(1- Af cost)) 2 


It should be noted that correlation coefficients a ljkl all 
have a Doppler-factor power of 4 in the denominator, 
which when multiplied by (Q/a^) 4 makes a ijkl 

proportional to the factctor k 4 . Shielding coefficients 
P ij depend on the number of turning points of g~(r,Q) 
as well as location of source point r o with respect to that of 
the turning point r a . For example, when there is only one 
turning point and r o <r a , we have 

Pu ~ exp{-2 1 -\/l g 2 (r)\dr } . 

Appendix B 

For a homogeneous isotropic turbulence correlation 
coefficients are related to the axial components with either 
Gaussian or exponential spatial functions 

/(£> = ex p(-7T<f // 2 ) and / (!) =exp(-7T<^/ 1) : 


I 

= I 

= I 


1 2222 

L 3333 

-Mill > 


■'1122 

= M133 

= 1 2233 

1 

> 

M212 

— ^1313 

~ I 2323 ~ 

Lj 

u 1111 ’ 


These relations hold for compact as well as non-compact 
source models. 
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Fig. 1 Energy spectrum for the exponential function /(£) = exp(-^ / f.) . 



0 0.05 0.1 0 0.05 o.i 


dx 2 (m) dx 2 (m) 

Fig. 2 Measurements of a two-point correlation (a) R 11 (^ 1 ), ( b ) R 22 (^ 1 ), (c) R n (£ 2 ), ( d ) R 22 {^ 2 ) 
in a high-subsonic jet. Eq. (2.14) with Exponential function (solid line); Gaussian (dashed- line). 
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Fig. 4 MGBK Noise spectra for Mach 0.5 cold jet. 

Predictions with model (a); compact source (solid 
line); non-compact source (dashed-line); data 
(symbol). 


Fig. 5 Effect of source non-compactness on spectral 

shape function of model (a). Compact 
source (solid line); non-compact source 
(dashed-line), (a) U / = 0.50, 1.0; 

(b) U / = 0.50 ; (c) V / = 1.0. 
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0.0001 




Fig. 6 Spectral shape function for the non-coinpact source 
with 

/(£) = exp {-K^ 1 / f) (solid line); 

/(|) = exp (-71% / (.) (dashed-line); 

(a) U I a x = 0.50 , (b) U / a x = 1.0. 



Fig. 7 MGBK noise spectra for Mach 0.5 cold jet with a non- 
compact source; 

Model (a) (solid line); Model (b) (dashed-line); 
data (symbol). 


Fig. 8a Comparison of the high frequency (HF) Green’s 
function with the adjoint Green’s function for a 
stationary ring source at indicated radial locations for 
a Mach 0.9 cold jet at St = 5.0. 
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Fig. 8b Green’s function comparison at St = 2.0. 


® 0.0015 
O 

O 0.001 

s 

^ 0.0005 


0.0015 

O 

O 0.001 


® 0.0015 

O 

O 0.001 




: r/D = 0.31 


Balsa, HF 
Adjoint 


40 60 80 100 120 1 



20 40 60 80 100 120 1 

r/D = 0.90 


20 40 60 80 100 120 1 

0 (from Axis) 


Fig. 8c Green’s function comparison at St = 0.90. 
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Fig. 9b Comparison of the axial cross-correlation coefficients 
of Tam et al. and MGBK. (i) Tam’s model (Eq. 3.11); 
(ii) MGBK model (Eq. 3.12) with (T = 0. 
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Fig. 9a Correlation coefficient I ]]u ( ) : MGBK (solid-line); Tam et al. (dashed-line). 


Fig. 9c MGBK cross correlation function (Eq. 3.12) 
with <7 = 0.40. 








Fig. 10 Spectral Shape function at 90°; MGBK (solid line); Tam et al. (dashed-line). 

(a) U / a„ = 0.325 , (b) U / a. = 0.65 . 



Fig. 11 Predictions for Mach 0.5 cold jet using a spectral shape function consistent with 
Tam et al. at 90°. 
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